Spin oscillations of the normal polarized Fermi gas at Unitarity 
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Using density functional theory in a time dependent approach we determine the frequencies of the 
compressional modes of the normal phase of a Fermi gas at unitarity as a function of its polarization. 
Our energy functional accounts for the typical elastic deformations exhibited by Landau theory of 
Fermi liquids. The comparison with the available experiments is biased by important collisional 
effects affecting both the in phase and the out of phase oscillations even at the lowest temperatures. 
New experiments in the collisionless regime would provide a crucial test of the applicability of 
Landau theory to the dynamics of these strongly interacting normal Fermi gases. 
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It is nowadays well accepted that a spin-polarized sam- 
ple of ultracold Fermi gases at unitarity, where the scat- 
tering length becomes infinite, is not superfluid, but 
rather a normal Fermi liquid provided the polarization 
is large enough (see, e.g.,.L] and the very recent experi- 
ment [5j and reference therein). The transition between 
the superfluid and normal phase occurs at the value 
{ni/nf)cc ~ 0.45 of the relative concentration between 
the minority- ! and the majority-^ component. Such a 
value is also known as the Chandrasekhar-Clogston limit. 

In an harmonically trapped geometry the critical value 
(ni/rif)cc is reached in the center of the trap at the 
value P = (JV t - N±)/(Nf + N i) - °- 77 of the (global) 
polarization of the gas, where Nf and iVj, are the paricle 
numbers of the majority and of the minority component, 
respectively. For smaller values of P we have a superfluid 
core in the center of the trap, surrounded by a normal 
shell. In the superfluid phase the densities are equal, at 
the border they present a jump, just beyond the border 
they have a ratio equal to the Chandrasekhar-Clogston 
limit and such a ratio decreases going outward in the 
trap. For larger value of the polarization no superfluid 
is instead present in the trap. The occurrence of this 
behavior is well confirmed experimentally [2|-l5j and the 
comparison with the theory of the normal phase, first 
introduced in [6j, is remarkably accurate concerning the 
critical value of P and the density profiles [1|. 

The theory developed in [l|, |6| is based on the assump- 
tion that the normal phase can be thought as a Fermi 
gas of polarons. A polaron in this context is just a 
single ^-impurity immersed in an ideal Fermi sea of J [- 
atoms, which behaves as a quasi-particle with an effective 
mass m* and chemical potential conventionally written 
as [i± = —3/5Afj,^, with A > since the interaction be- 
tween the two species is attractive. The values of the 
dimensionless parameters A and m* jm have been calcu- 
lated by means of different theoretical approaches |6T-fl"ol 
and first measurements have been reported in [3j-l5j, lll| . 

In the presence of an harmonic trapping V/i D (x) = 
l/2m{uj 2 x 2 + uj 2 y 2 +ulz 2 ), the polaron is well described 
by the single-particle effective Hamiltonian Qj 




showing that the impurity feels the renormalized oscilla- 
tor frequencies 



(2) 



Accordingly it will oscillate with frequencies renormal- 
ized by same factor \/(l + 3/ 5 A) m/m*. When more 
impurities are present their dynamic behaviour will be 
affected in a different way by the interaction with the 
majority component and the frequency of the oscillations 
will depend on the value of the polarization. 

Very recently the results of such an experiment have 
been reported in [3j, where the frequencies of the ax- 
ial breathing modes for the majority and the minority 
components have been measured at different values of P. 
Intuitively one expects that there are both in phase and 
out of phase modes of the two spin components, the lat- 
ter corresponding, in the limit of zero concentration, to 
the polaron oscillation. Indeed in the experiment [3j it is 
seen that the majority and the minority component share 
a common frequency (in phase mode), but the minority 
one presents also a second branch at higher frequency. 

The aim of the present work is to study the collec- 
tive mode frequencies of the unitary normal phase as a 
function of the polarization. To this purpose we develop 
the Landau formalism of quasi-particles to describe the 
collisionless regime of the polarized Fermi gas. 

A first important issue is the problem of the damping of 
the oscillation, since collisions can be very effective in this 
strongly interacting system. Let us consider highly po- 
larized Fermi gases at a temperature T smaller than the 
degeneracy temperatures T^ a of both the components. 
The damping of the dipole oscillation was calculated in 
Il2] and can be written as 
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where u>q is the oscillation frequency of the order of the 
harmonic trap frequency, SX is the amplitude of the os- 
cillation and i?-j- the Thomas- Fermi radius of the majority 
component. The calculation has been carried out for an 



homogeneous system and then applied to a trapped gas 
via the local density approximation. Aside from some 
numerical factors expression ([3]) is typical for the life- 
time of quasi-particle in a Fermi liquid [13| |. As can be 
inferred from ([3]) the conditions for reaching the collision- 
less regime loqt 3> 1 impose severe conditions from the 
experimental point of view. Indeed in the experiment [3| , 
where the axial compressional mode was investigated in 
a highly elongated trap, the oscillation of the minority 
component exhibited strong damping and its frequency 
could be extracted only through the Fourier transform 
of the signal. Moreover also the in phase axial breathing 
mode, which is dominated by the majority component, 
was found to approach the collisionless value of the fre- 
quency only for very high polarization [3j . 

The collective modes in the collisionless regime can be 
derived by using the variational principle SS = applied 
to the action integral 



(48-/V-J-) 1 ' 6 ^Jft/jj /muji, i = x, y, z with w 3 = uj x u> y uj z . In- 
stead the radius of the minority component is quenched 
with respect to the non interacting gas due to the at- 
tractive nature of the force. By taking a Thomas-Fermi 
description for the minority component (which holds with 
good accuracy for a large class of experimentally available 
configurations), one finds that the radius of the minority 
component is given by the simple expression 



S = I dt(V\H - ihd t \^) 

= f dt(E - <tf|»m|tf», 



(4) 



where | Vt) is a multi-parameter many-body wave-function 
and E = (&\H \i$f) is the energy functional of the system, 
that we write in the form 
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where T a j1m is the kinetic energy density of the species 
a. The functional (JSJ accounts for the interaction be- 
tween the majority and the minority component through 
the local term in A and the last integral, which is nec- 
essary in order to keep into account that the polarons 
(4,-particles) acquires an effective mass due to interac- 
tion effects and that galilean invariance implies the pres- 
ence of a counter current term jl of the majority com- 
ponent. We include the latter term for completeness, 
although, since the effective mass is just 10 — 20% larger 
than the bare one and the effect of the current term on 
the frequencies turns out to be fairly small. Expression 
([5]) corresponds to a typical energy functional to be used 
in time-dependent Hartree-Fock approaches in the con- 
text of small amplitude and low energy oscillations. This 
approach is equivalent to Landau's theory of Fermi liq- 
uid. 

At equilibrium the kinetic energy density in Eq. ([5]) 
reduces to r a = f?{^-K 1 n G ) 1 ^ n <y an d j-f = and, thus, 
the energy functional Eq. §5$ can be used to calcu- 
late the density profiles using standard variational pro- 
cedures. The calculation at equilibrium shows that at 
very high polarization the majority component is scarcely 
affected by the interaction and, in paricular, its ra- 
dius is in practice given by the ideal gas value i?2 i = 
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From Eq. ([6]) it is seen that the minority radius is quite 
flat as a function of P except at very high polarization 
when it goes to zero. We find that this behaviour is 
reflected in the collective mode frequencies (see below). 
In order to determine the solution of the time- 
dependent problem we develop a variational approach- 
based on the scaling transformation 



ip<r(r,z,t) 
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applied to the single particle wave-functions t/v of the 



two spin species a =~[, \,, and where r 



y and z 



are the radial and the axial coordinate, respectively. The 
scaling transformation depends on 4 + 4 time-dependent 
parameters and the corresponding equations are obtained 
by imposing the variation SS — with S as in Q with 
(^|if|^ r ) given by Eq. (J3J) - With respect to a typical 
hydrodynamic energy functional, Eq. ([5]) accounts for 
the deformation of the Fermi surface produced by the 
scaling ansatz. This effect arises from the kinetic energy 
density term T a and exploits the elastic nature exhibited 
by a Fermi liquid in the collisionless regime. The use of 
hydrodynamic theory would actually yield wrong predic- 
tions for the oscillations of these Fermi gases. 

By making the choice Q for the scaling transforma- 
tion we implicitely assume that the equilibrium configu- 
ration is axially symmetric and restrict our study to the 
coupled surface compressional modes of axial/radial na- 
ture. The reason why we restrict the analysis to such 
modes is related to the fact that the easiest experimental 
way of exciting the spin modes in trapped Fermi gases is 
through a sudden change of the value of scattering length 
(see Appendix A for the single polaron case). This proce- 
dure, which mainly affects the motion of the impurities, 
is not able to excite other important oscillations like, e.g., 
the dipole mode of the minority component and it is the 
actual procedure employed in the recent experiment re- 
ported in Ref. [3]. 

The collective modes are just small oscillations around 
equilibrium, i.e., solutions of the equations of motion de- 
rived from the action expanded to second order in the 
scaling parameters. The first order expansion of the ac- 
tion S with respect to the scaling parameters takes con- 
tribution only from the energy functional Eq. ([SJ . Then, 
the condition SS = provides a relation between the ki- 
netic and the interaction energy equivalent to the virial 



theorem. Defining the effective mass as m/m* = (1 + a) and the averages N a (f) a = J f(r, z)n a (r, z), we obtain 
I 
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Similar relations hold along z. In Eq. © all the densities 
are to be calculated at equilibrium. The term in the 
action which depends on the time derivative of the wave- 
function does not give rise to linear terms due to time 
reversal symmetry. The quadratic term is given by 

(V\ihd t \*)W = 2Y / Na((r 2 ) a a^ a + (z 2 ) a /3 a x*)- (10) 



Summing up all the contributions and imposing the vari- 
ational procedure 5S = 0, we get eight coupled equations 
of motion. Four of them represent continuity equations 
and are relations between the current parameters £ CT , Xo- 
and the densiy ones a ff ,/3 ff . In this way we are left to 
solve a linear system of 4 equations. 

We find that the two lowest frequency are almost in- 
dependent of the ratio N^/Nf and they are very close 
to the ideal gas values ui = 2u± and u = 2uj z . We 
call them the in phase modes since the majority and the 
minoity components move in phase. The frequencies of 
the other two modes, that we name out of phase or spin 
modes, can be written as u — 2CiUJ± and u> — 2C*2W Z 
with the renormalization factors turning out to be very 
close, i.e., C\ ~ C2. They basically correspond to the 
radial and axial motion of the minority component mov- 
ing in opposite phase with respect to the majority one. 
In the limit of a single impurity we recover the values 
d = C 2 = ^ (1 + 3/5 A) m/m* of Eq. ©. 

In Fig. (TTJ we report the result for the axial spin mode 
as a function of the polarization of the system together 
with the experimental data of [3j. We notice that at 
high polarization the correction to the polaron frequency 
Eq. ([2]) , as a function of the polarization, follows the law 
(1 — P) 1 ' 6 characterizing the radius of the minority com- 
ponent (see Eq. ©). From Fig. (TTJ) we see that there is a 
qualitative difference between the present theory and the 
experiment Q, but a remark is due here. Our prediction 
was derived in the colliosionless regime, while the exper- 
iment clearly shows that the in phase mode frequency is 
strongly affected by collisions even for the highest polar- 
ization available. Thus the question of how the observed 
frequency in Fig. ([T]) can be compared with our pre- 




diction has not an obvious answer [14J. Our prediction 
agrees better with the experimental datas at the high- 
est polarization points, where the collisionless approxi- 
mation is better satisfied [3[. In this respect we expect 
that a measurement of the radial compressional mode, 



FIG. 1: Frequency of the axial compressional mode as a func- 
tion of the polarization P. Dashed line: the single polaron 
mode frequency. Solid line: the collisionless value of the mode 
frequency obtained via the variational principle described in 
the present work. Points: experimental data reported in [3|. 



should give a much better insight into the problem since 
in highly elongated traps the radial frequency is much 
higher than the axial one and thus collisions are less ef- 
fective. Such an experiment would provide a stringent 
test for the applicability of Landau theory to the dynam- 
ics of these novel systems. 

As a last comment, we should remind that the energy 
functional ([5]) is strictly valid only at high polarization, 
being an expansion in the concentration n^/n^ [l|,|6J|. It 
works however quite well when applied to investigate the 
static properties even till the critical concentration. For 
instance the inclusion of the next term proportional to 
(n^/n^) 2 [l|, [Hj, [l5[ introduces only small corrections to 
the values of the radii, even close to the critical polariza- 
tion limit Pc — 0.77. Such a term does not change the 
picture emerging from the present analysis. 

A natural extension of this work would be to study the 
collective mode frequencies for the impurity out of reso- 
nance. In particolar for positive values of the scattering 
length it would be interesting to study the distinction 
between the polaron state and the molecular one from a 
dynamical point of view. 
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Appendix A: Sudden change of the scattering length 
and polaron dynamics 

What does it happen to the dyamics of a single impu- 
rity in a Fermi sea if one suddenly change the value of 
the scattering lenght? 

In presence of an harmonic trapping potential Vho we 
can write for the polaron an effective time dependent 
Hamiltonian as 
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where m(t) is the istantaneous value of the effective mass 
and we define an effective frequency u)f(t) — m/m(t)(l + 
3/5A(t))w|. The Hamiltonian (|A1|) . is easily diagonal- 
ized to the form H = J2 hu>i(t)ai(ty Oi(t) by introducing 
the istantaneous ladder operators 
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where Ui (t) = y/m(t)u>i(t) is fixed by the way the scatter- 
ing length is changed in time. In the Heisenberg picture 



ai(t) = -iui(t)a,i(t) + -^-rxaiitf ', 

Ui{t) 
Ui{t) 
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The first terms on the RHS describes the trivial evolution 
of the operators, while the second ones take into account 
the explicit time dependence of the parameters in the 
Hamiltonian and are responsible for the mixing between 
di and a], which eventually depends on the rate at which 
the scattering length is varied in time. 

Once the functional form of Ui(t) is fixed it is easy 
to calculate the evolution of the width of the impurity 
wavefunction Axf. In the case of a sudden change (see 
also [16[) it is easy to write the evolution of the ladder 
operators as 
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(A6) 
where u>j/ is the final value of the effective trapping fre- 
quency, Uifi and Uij are the initial and the final val- 
ues of Ui, respectively. Let us suppose we start from 
a situation where we do not have any interaction be- 
tween the two-hyperfine levels for which m(0 _ ) = m and 
A(0~) = 0, for which Ui t o — ^/muji, and the system is 
in its ground state. Suddenly we change the scattering 
to a very large value (unitarity), for which m(0 + ) = m* 
and A(0+) = A, i.e., u hf = ^UTi(mm*{l + 3/5A)) 1 / 4 . In 
this case the width of the polaron along the i-th direc- 
tion will oscillate at a frequency given by twice the oscil- 
lator frequency ^ around the new equilibrium position 
Axf = l/(4mwi) (l + (m*/m(l + 3/bA)^ 1 ) according to 
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Thus, the measurement of the frequency of the oscillation parameters of the polaron. 
and of its amplitude provides the values A and m* of the 
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